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ABSTRACT 


The application of error detection and correction codes 
has advanced markedly with the advent of digital technology. 
However, the strides made towards employing the techniques, 
encoding and decoding, properly have been rather limited 
by the Naval forces in the United States. This paper deve- 
lops a computer program, which if utilized properly, would 
aid in deciding what error correcting scheme is best suited 
for a specific channel. 

The results obtained from testing a rate 1/n convolu- 
tional code, over a simulated channel, using a Viterbi 
decoder shows that this is an effective analysis procedure. 


Though the test runs were lengthy, much of the time required 


was for noise simulation. This would not be a factor if actual 


channel noise recordings had been available. 
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I. INTRODUCTION 


The distortion of a digital transmission, caused by 
noise, may require that techniques be used allowing for the 
detection and eventual correction of errors. In digital 
communications, the most significant performance parameter 
from either a bit or message standpoint is that of probabi- 
lity of error. A means of minimizing the probability of 
error is therefore essential to effective communication. 

Error detection and correction was given great impetus 
by Shannon's paper of 1948, [Ref.7 ], which extended the 
promise of reliable recovery of digital data perturbed by 
noise (Shannon's Second Theorem). The noisy coding theorem 
provides that messages can be transmitted with arbitrarily 
small error if the source rate is no more than the channel 
capacity.! While this sets a goal and the conditions neces- 
sary to attain it, the precise method for obtaining a speci- 
fic scheme was not set forth. However, Shannon did use a 
random coding scheme in proving this theorem. As was shown 
by Rice in 1950, [Ref. 6 ], following Shannon's example, 
choosing codewords randomly leads to the result that as the 


code becomes very long, channel capacity is approached and 


l Channel capacity is the maximum average information 
input to a channel, which means it is properly matched to 
the channel less the average uncertainty at the receiver 
due to channel noise. 
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the probability of error can be made to decrease exponen- 
tially with code length. But choosing codewords randomly 
is not a practical scheme. 

Until the 1960's the application of channel coding theory 
was slow in its development. The establishment of digital 
circuit technology and the realization of a theoretical 
channel (satellite communication; AWGN) provided a powerful 
stimulus for the utilization of practical error detection and 
correction techniques. The technological improvements have 
continued until today fairly large memories are held by 
small chips, and the components needed to implement powerful 
coding schemes are available in reduced sizes and at reason- 
able prices. 

The objective of this paper is to demonstrate the bene- 
fits derived from a minicomputer (DEC PDP-11/40) simulation 
of channel noise applied to a specific error detection and 
correction scheme. This simulation allows for the determina- 
tion of a probability of bit error for a specifiable noise 
density. Results from the simulation are then used to deter- 
mine statistical trends of a particular error correction 
technique for various internal coding parameters. 

The error detection and correction scheme chosen by the 
author is the Viterbi decoding, [Ref.3 ], algorithm on a rate 
1/n convolutional code. Both the encoding and decoding 
technique are of current interest in extraterrestrial com- 
munications, with the Viterbi algorithm being principally 


considered for satellite communications. Therefore, in 


11 





order to provide a complete presentation, the scope of this 
thesis will encompass the basic principles of convolutional 
encoding and Viterbi decoding, as well as the structure of 
the simulation program, results and conclusions. 

All programming was accomplished on the DEC PDP-11/40 
in machine language, thus presenting the opportunity to 
gain valuable insight into the computer structure while 


completing this thesis. 
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II. BACKGROUND 


Digital communication systems (figure 1) are usually 
designed to minimize the probability of error of the re- 
ceived data (bits), introduced by noise (gaussian, burst, 
fading, multipath, etc.). Thus, error detection and correc- 
tion may be varied according to the degree and type of noise. 
There are two common forms of error protection: (1) retrans- 
mission and (2) controlled redundancy. The first may be 
applied to systems which are not critically affected by the 
ensuing delay, whereas the second does not require signifi- 
cant delays. This thesis deals with the latter in an 
attempt to gain statistical knowlodge of the particular 
coding scheme (controlled redundancy) under study. 

Controlled redundancy techniques are commonly divided 
into two groups: (1) block codes and (2) convolutional 
(tree) codes. Again, owing to the characteristics of the 
channel, one of these two codes may be chosen, along with 
a suitable decoding scheme, to achieve the desired probabi- 
lity of error. A decoding scheme for block codes is batch 
oriented, as is its encoding, using definite algebraic 
operations related to the segmented structure of the code. 
On the other hand, convolutional codes are decoded by a 
statistical procedure due to its continuous (bit by bit) 


nature. 
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A. PRINCIPLES OF CONVOLUTIONAL ENCODING 

In 1955, P. Elias first proposed the use of convolutional 
(tree) codes for the discrete time memoryless channel, 
[Ref 2. This technique extended the promise of providing a 
class of codes (linear) whose performance would prove supe- 
rior to that of block codes of the same length, [Ref.9 ]. 
The development of these codes also gave promise of provid- 
ing for a decoder complexity increasing no more than linearly 
with block length and/or encoder memory. These conjectures 
have been verified, for the most part, by such contributors 
as Fano, Reiffen, Forney and Berlekamp. [Ref. 4] 

Since convolutional codes form a definite discipline 
in coding theory, there are some aspects that will require 
clarification in conjunction with the implementation of the 
code. The next few subsections will accomplish this task. 

T. Description 

The process of block encoding segments (blocks) of 

an indefinite length input sequence into code blocks pro- 
ceeds in a discrete manner, i.e., the contents of one block 
has no bearing on the encoding of another block. The code 
blocks may be completely changed from the input block or 
they may include the input block with check bits determined 
by a definite algorithm. However, when an input sequence 
is processed through a convolutional encoder the input/ 
255 Continuous (sequential bit to bit dependence in 
the encoder) and the output is generated at a specifiable 


number of code bits for each encoder input bit. 
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The term "convolutional" originates from the observation 
that the encoded sequence can be regarded as a convolution? 
of the input sequence with related generator sequences. An 
example of such a procedure is the binary convolution 
(denoted by *) of two sequences, x and g, where (+) and ° are 
binary addition and multiplication, respectively. 


Imput sequence: 


x (tg) = (xgXiX4X4X, cv... UĞUR 


Generator sequence: 
R g (tg) = (9998195959, .. — = (1 0 0 l l o 5. , 
where all bits in x and g are zero for x. and 9... 


Then the output (code) bit formed by the convolution of 


x (t9) with g (to? is 
y (to) = x(tg)* g(tg) =0. 


For time tg t At, the sequences and output are 


x(totAt) = (xgxixj4X4X, 277 — Clo le O ne), 
g(tgt^t) = (9091929394 mu C130 0 ds lis 
y(t tat) = (xo . go) = (1 ` 1) = 1. 


For time to + 2At, the sequences and output are 


?The convolution of two functions, fi and fo, is equal 
to 


| e, (t-1) f, (ar, 


which means one function is folded in time, then their 
product is integrated as one is shifted pass the other. 
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x(tgt2At) = (x9X1X4X4X, 235 “ (L O 1170”), 
g(t,¿+24t) = (9091929394 29 (1 0 0 1 1 ...), 
y(tot2At) = (Xp . O, (xi ° 90) = 0. 


And the same procedure is applied to subsequent increments 
0:77: (one bit shift). 

This technique extends the concept of block encoding 
to permit memory from bit to bit continuously as against 
memory within the block. Using the operation shown above, 
a conventional convolutional encoder (figure 2) may be 
defined as a linear sequential machine with k-inputs and 
n-outputs, where n » k, usually. This machine is constant, 
linear, causal, and finite state with operations over a 
finite field F, commonly binary. 

The parameter k is the number of input bits which 
enter the encoder in a time increment (At). These bits 
along with others retained in the encoder memory (finite 
length), form n-output bits (code bits). For the remainder 
of this paper k will have a value of 1, but it should 
be remembered that this parameter can take on other values, 
as can be seen in figure 3. 

The definition specified a finite state machine, where- 
as the example dealt with indefinite length sequences. 
Since an infinite length memory is not practical it becomes 
necessary to limit the basic concept (generator sequence) 
and define some new expressions employed in the actual 


encoder configuration. 
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Figure 2. Conventional (n,k) Convolutional Encoder 
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a. Code Rate (k/n) 

The code rate is an expression of the number 
of k-input (1) bits per n-output bits, that is, k/n (1/n). 
This value may be considered a very close approximation 
to the actual value. The actual value would also include 
a very small number of zero code bits used to terminate the 
code sequence. 

b. Memory Length (M) 

For a practical encoder g must be limited in 
length. The oc length is a measure of this length and 
is the minimum number of memory cells required to generate 
a code. A representation of an encoder, with a memory 
length of 2, is shown in figure 4. This configuration 
corresponds to the Mealy machine in automata theory. 

c. Constraint Length (K) 

A more conventional expression of the length 
of time which an input bit affects the output sequence is 
that of constraint length. “This term is simply the memory 
length plus one, and is depicted as the number of memory 
elements in a Moore machine configuration of an encoder 
(figure 5). As may be expected this value is one factor 
in determining the complexity of the error detection and 
correction scheme. 

d. Free Distance (de) 

The term distance, as applied to coding, refers 

to the number of differing bits between two code sequences 


of the same length, as seen below: 
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code sequence 1: (1 00 10 11 00 1) 

code sequence 2: (101100100 1) 
Differences between the two sequences exist in the third 
and sixth bits only, therefore the distance between the two 
sequences is 2. When using block codes to construct a code 
sequence the minimum distance between all pairs of codewords 
(code blocks) is a definite indication of the error detect- 
ing and correcting capability of the code. However, if the 
code 7 not employ a block configuration, the use of this 
value may not be justified to indicate the codes error 
detecting and correcting capabilities. 

This is the case when discussing convolutional codes. 
Free distance is the term for the value expressing the 
theoretical error correcting capability of a convolutional 
code. This value is defined as that minimum number of set 
bits occurring in a code sequence which resulted from the 
input of a nonzero sequence. Table 1 shows how the free 
distance (дӱ) сап be determined for a rate 1/2 convolu- 
tional code with K=3. As can be seen, the larger the 
value of K then the number of possible input sequence to 
be considered also increases at a great rate. But with the 
aid of computers the tedious grinding process is accomplished 
rather quickly. 

The convolutional encoder is constructed to take full 
advantage of the free distance applied to a specific code 
rate and constraint length. This is accomplished through 


the proper selection of generator sequences. However, it 
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should be noted that with the computers currently available 
the drudgery can be passed on to the machine. 
2. Properties 

The properties of the encoder, as related to the 
generator sequence, are discussed in this section. For 
convenience, the set of generator sequences äi e Bal 
used to generate the convolutional code is denoted by 
[G], a matrix representation of the encoder discussed in 
the next section. Figure 6 is useful in relating the 
following discussion to a physical communication system. 

a. Property 1 

Foremost, the useful encoder should generate 

a code which will yield the fewest number of errors in the 
codeword estimator. The estimator is a demodulation scheme 
at the receiver determining as accurately as possible the 
received sequence before actual decoding takes place. Along 
with this error minimization, the complexity of the estima- 
tor requires minimization to develop a useful error correc- 
tion system. Error minimization and complexity minimization 
together may require some compromise in the practical 
system depending on the desired reception quality and sys- 
tem environment. 

b. Property 2 

The encoder [G] must have an inverse relation 


for decoding purposes. This is realized by the relation 


1 


yid == xtejtc) L = pP 


bs 
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=L 


for all x. The matrix (decoder), [G] is physically 


realizable and is a pseudo inverse, in that, the expression 
E –1 : : 

[G] [G] ^, yields a delay DP, where p is the number 

of time intervals before decoding occurs in the expression. 


If p is zero the decoding occurs on receipt of the 


sequence and the above relation becomes 


y[G] l - x(G)1G17İ = x, 


EJ 1 


where [G] SI лат 1 — 0), 
A mention of the possibility of catastrophic error 
propagation, [Ref ], is required at this time to stress 


zü 


the point that te) must be feedback-free. This may be 


interpreted as meaning that the n generator sequences, 
in polynomial form should not have a common factor. If a 
common factor exists then the input of a sequence with a 
finite number of set bits to the encoder, may be decoded, 
after noise is added, as a Sequence with an.infinite 
number of set bits. This is catastrophic error propagation 
and is avoided by making [G], thus Ton feedback-free. 
c. Property 3 

The matrix [G] must also meet the requirements 

mentioned in the definition of a convolutional encoder. 


(1) Constant (time invariant). The time 


invariance of the encoder is represented by 


[G] (pPx) = pPy , 


2? 





which means if all the inputs are shifted in time, then 
all the outputs are shifted accordingly. 

(2) Linear. The output sequence resulting 
from the superposition of two input sequences, x, must 
be equal to the superposition of the two output (code) 
sequences, y, that would result from the inputs entered 
separately. This is also necessary for the multiplication 
of an input by a scalar as shown below. If G:x>y, then 
G(x, +x.) = G (x4) t G(x.) = y} + Y), and G(ax;,) = aG (x, ) - 
ay, where o is an element in the field F (usually GF(2)). 

(3) Causal. The existence of a nonzero output, 
y, prior to the input of a nonzero Sequence, x, is for- 
bidden. This is accomplished by the encoder output being 
zero when its memory elements are all zero, therefore zero 
inputs. 

(4) Finite State. The states of the encoder 
are finite in number due to the value of the constraint 
length. Thus, a binary encoder has 2. possible states. 
Each state being those input elements involved in the 
generation of the next output (code) bit, vhen combined 
vith the next input bit entering the encoder. 

d. Property 4 

The chosen encoder (G) should have the minimum 

number of memory elements to generate the code. 
e. Property 5 
The chosen encoder İGİ should generate a code 


yielding the fevest number of decoding errors per error 


¿8 





event (channel noise).  Thereby, the fastest correct 
decoding decision is made at the receiver. 

The preceding paragraphs provide a basis by which the 
acceptable class of encoders is greatly reduced in number. 
However, the selection of [G] in its optimum form for a 
given channel is a long drawn out process, possibly best 
Suited for computer analysis. Theoretical selections have 
been made by many people, [Ref. |, who based their deci- 
sions on the free distance alone. 

3. Representation 

The encoder configurations used to this point are 
accurate physical descriptions, but they lack a convenient 
form needed for analysis. There are five encoder repre- 
sentations which show varying degrees of the code's struc- 
ture. Each is discussed in the following sections. 

a. Polynomial 

Using the D-transform, the input and output 
sequences appear as polynomials of degree much greater then 


the constraint length (K). 


É -ı 2 
mput, X(D) = ... + x ,D + Xo + x,D + xD + 


z i 
Output, Y; (D) mo... + e + Y; O + y; jP 


2 
+ Yo + ..os 


where 1 s i < n. The generator sequence, for the Mealy 
machine configuration in figure 7, is a polynomial of 
degree equal to K-l. The mapping of X + Y. is accom- 


plished by polynomial multiplication. 
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where L is the bit length of the input sequence. This 
representation has little advantage in a code structure 
analysis, but it is very definitive as to the procedure for 
code generation. 
b. Matrix 

A matrix, [G] representation of a rate 1/n 
encoder starting at state zero and time zero is achieved 
by the manipulation of each generator sequence. When 
y =[G]x ,[G] may be shown as a matrix whose row vectors 
are the generator sequences (figure 8a). Using this nota- 
tion restricts x to be a length K. However, the require- 
ment to change x for each code computation is eliminated 
by the matrix in figure 8b. This notation incorporates 
the shifting of x past each g.(1 < i < n) into the 


matrix  [G]. For example, (figure 8 follows the example), 
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This form allows for a fast generation of the code and 
is definitely a model which could be implemented on a com- 
puter with little difficulty. Again, there is a lack of 
BBB Һе inte the code structure except for the possible 
application of matrix algebra theory. 
cı Tree 

The most common representation is a tree dia- 
gram (figure 9) vhich incorporates the branch and nodal 
properties of the code generation. The base branch corres- 
ponds to the initial state of the encoder prior to a non- 
zero entry. The first node (a) is state zero for the en- 
coder. When a 1 begins the sequence the first lower 
branch is chosen and the code is found on this branch 
leading to node (b) or state 1. The subsequent input 


bits dictate whether the tree is followed up or down, 
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corresponding to inputs of 0 and 1 respectively. If 
this is continued to the Qum branch, all of the possible 
input sequences of K bits have been represented, thus 
each possible encoder state. Beyond the — branch it is 
readily recognized that the code symbols on the branches 
leaving from the tvo nodes labeled (a) are identical iden- 
tical to those at the base of the tree. Even though the 
structure has become recurrent, for long sequences this 
diagram may become too large to analyze easily. 

d. Trellis 

Since the input sequence (100 ...) and 

(000 ...) generate the same n-bit codewords after the 
third branch (figure 9), then both nodes labeled a can be 
joined together. The recursive character of the tree 
diagram lends itself to be redrawn with remerging branches, 
thus forming a trellis (criss-cross) diagram (figure 10). 
The encoder state is the basis for the trellis diagram. 
These states represent those input bits in the memory cells 
which will generate the code bits when the next input bit 
enters the encoder. Therefore, there are Se? states on 
either side of the diagram representing the present and 
next states. The branches connecting the states denote the 
code bits being generated when a 1 (dashed line) or a 
0 (solid line) enters the encoder. Again, after the third 


input bit the complete trellis diagram is specified and is 


repeated for all succeeding branches. 
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For short constraint lengths (K < 10), this repre- 
sentation is highly desirable due to its compactness. The 
trellis diagram is used in the Viterbi decoding algorithm 
and easily lends itself to be listed in a computer program. 

e. State Diagram 

The last encoder representation is similar to 
the trellis, but is spread out more to allow more parameters 
to be placed on the branches between states (nodes). The 
state diagram (figure 11) has the input 0 and 1 denoted 
in the same manner as for the trellis, but some additional 
notation is placed on the branches. 

(1) pi, The value of q denotes the weight 
of the code bits Eon that branch. 

(2) Ce The value of r denotes the branch 
length of a path from state x to state y. 

(3) N”. The value of s denotes the number 
of input one branches encountered in a path from state Xx 
to state y. 

All of the above notation is used in figure ll and the 
results of a state path for the input sequence (10 1 0 0) 
are shown. 

Since the state diagram is a directed graph, a transfer 
function can be determined using the theory of directed 
27 5:5 [Ref. 4, pp 2A-1-2A-11]. The transfer function con- 
sists of various powers of the three measurements listed 
above. These values are used to determine the properties 


(D, L, N) for all paths in the convolutional code, when the 


transfer function is in expanded form. 
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Figure 11. State Diagram Representation of 
Convolutional Encoder in Figure 5. 
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The trellis and state diagrams are the most complete 
representations and offer the analyst a chance to study 
the path structure of the convolutional code. Each diagram 
yield a compact representation, easily arranged for small 
values of K. However, the looping branches necessary for 
the construction of a state diagram become very confusing 
and difficult to arrange for larger constraint lengths. 
Therefore, the trellis diagram is not only easy to conste 
and understand, but if the values of D, L, N, are needed 
for analysis they can incorporated into the diagram in an 
orderly fashion. 

The discussion of convolutional encoding principles 
was basic and is by no means a complete detailed course in 
convolutional codes. However, the purpose of this section 
is to introduce the reader to concepts which will aid him 


in his understanding of the thesis computer program. 


B. PRINCIPLES OF VITERBI DECODING 

Decoding is the inverse operation to encoding and is 
intended to recover the source bits with all, or almost 
all, of the channel errors removed. The decoder, which 
may be implemented as hardware or software, utilizes the 
encoded bits to detect and/or correct errors. Error detec- 
tion is similar in complexity to the encoding operation. 
Error correction, however, is inevitably a more complicated 
process than encoding, since the goal is to reduce the 


probability of a decoded bit error. 
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A scheme for decoding convolutional codes was proposed 
by Viterbi in 1967, [Ref.8]. It was shown by Viterbi that 
an optimum decoding procedure existed for a statistically 
independent (bit to bit) input sequence transmitted over a 
channel whose errors occur independently from channel bit 
to channel bit. Subsequently, Forney found that the Viterbi 
decoding algorithm is synonomous with maximum likelihood 
sequence decoding, [Ref. 3]. 

1. Description 

Simply stated, the Viterbi algorithm is a solution 
to the problem of finding the most likely encoded sequence 
through a state (finite) diagram representing the encoder. 
For a decoder to minimize the overall error probability 
of a decoded bit by brute force, maximum likelihood decoding 
would mean calculating the likelihood of the received 
sequence on all paths of the encoder state diagram.  How- 
ever, there are two factors which reduce the complexity of 
this problem. The first is the fixed periodic structure 
of the encoder trellis diagram, and the second is the code 
characteristic of remerging paths after the same K input 
bits are applied to two different paths (figure 12). For 
these reasons, the trellis diagram is an ideal tabular 
representation of the flow of the code at any instant. 

The decoding process, as mentioned above, was found 
optimal for a area independent input sequence in 
discrete time. Along with this stipulation, the channel 


noise is also memoryless, as is the case for a binary 
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symmetric channel (BSC) with hard decision demodulation 

(21 levels) or an additive white gaussian noise (AWGN) 
Channel with soft decision demodulation (29 levels, Q > 1). 
The BSC errors transforma 0 toa 1 anda 1 toa 0 
and occur independently from bit to bit with probability 

p. In the AWGN channel, the probability of a given quan- 
tized value (0 to 29-1) Of a received bit is determined 
from the gaussian probability density function. These 
values of the received bits are used to determine the most 
likely received sequence. 

A scoring procedure is employed indicating the trellis 
path which shows the least difference from the received 
sequence. This is accomplished by using the concept of 
distance for the BSC and a numerical difference for the 
AWGN channel. The latter may be called the innerproduct 
of the received bit and the calculated (trellis) trans- 
mitted bit. This is implemented digitally by taking the 
exclusive OR of the two Q bit representation of the bits 
above. Both (BSC and AWGN) techniques yield a value that 
can be used to determine the most likely (lowest score) 
path. 

2. Implementation 

In order to describe effectively the implementation 
of the Viterbi decoding algorithm an example follows with 
step by step explanations supplemented by appropriate dia- 
grams. This example will be for a rate 1/2 convolutional 


code with K = 3. The generator sequences are 5g Q 0 1) 
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and 79 (l1 1). The following sequence = the encoder in- 
Ell 0 1107000 1 1011000 0 10...). 
The corresponding code sequence is 
2177) 00710 10 11 00 11 10 10 00 01 11 00 11 01 11 00 ...). 
An error sequence denoting BSC noise perturbs this sequence 
prior to being received at the decoder. The error sequence 
is (01 00 00 00 01 00 00 00 00 10 00 00 00 00 00 10 00 00 
...). Then, the received sequence is (10 01 00 10 11 11 00 
11510 00 00 1] OO 11 11811 00... 8 
The following subsections are the steps of the Viterbi 
algorithm decoding the perturbed (received) sequence. 
a. Step 1 
The following figure and tables are set up to 
indicate the trellis structure, the path (survivor) sequences 
and the path scores. 
(1) Trellis. (Figure 10). 
(2) SSEQ(I). The present state survivor 
sequence. 
(3) SCORE(I). The present state survivor 


sequence score. 


(4) SSEO(J). The next state survivor sequence. 
(5) SCORE(J). The next state survivor sequence 
score. | 
b. Step 2 j 
An initialization of tables (2) and (3) is 
shown in figure 13. Tables (4) and (5) are not considered 


until the first group of 2(n) received bits enter the 
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decoder. The all zero sequences placed in SSEQ(I) depicts 
the assumption that the encoder had no set bits entered 
prior to the beginning of the message under consideration. 
A decoder parameter is the length of the survivor sequences 
used to hold path estimates until a bit is decoded with the 
desired probability of error. This length is called the 
decoder constraint length (DCL), whose value is 6(2 x K) for 
this example. Also initialized is SCORE(I) so that a score 
of Zero is in state 0 to denote that the encoder is 
assumed to have started from this state. The other states 
are assigned scores which demonstrate the unlikelihood of 
the encoder starting from these. Normally, a value of 
2 x K is sufficient, which is 6 in this example. 

c. Step 3 (Figure 14) 

This step is the first in the actual decoding 
process. However, lt will become obvious that it is recur- 
sive throughout the remainder of the alqorithm. Now, the 
first 2 received bits (10) enter the decoder. These bits 
are compared with each pair of code bits on the branches 
of the trellis and a distance (A) is determined for each 
branch. This A is then added to the present score of the 
corresponding present state. Now there are two scores 
coming in on the branches to each next state. At this time, 
a decision is made as to which branch, coming into each next 
state, has the lowest score. If the scores are identical 
then an arbitrary choice can be made, decode estimate is 


zero for this example or an alternative may be exercised to 
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Present State Decoded Bit = O From State 0 


Figure 14. Next State Tables at the Completion of the 
First Shift of 2-Bits Into the Decoder. 
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further evaluate the A added to the score in SCORE(I). 
When all the next states have a score, they are entered in 
SCORE(J). Then the next state's survivor sequences are 
determined for entry into SSEQ(J). This is accomplished by 
determining what bit, 0 or lwas dropped from the present 
state at the end of the minimum score branch for each next 
state. Then, this bit is shifted into the beginning of the 
appropriate present state sequence (SSEQ(I)) and the entire 
sequence is then placed in the SSEQ(J) table, corresponding 
to its next state. 

The oldest bit just shifted out of the sequence in 
SSEQ(I) at state 0 is the most likely decoded bit since 
it is associated with the minimum A (0) sequence in the 
present state. The bit decoding, after a length of DCL, 
minimizes the effect of randomly spread errors in the 
received sequence. Again, if the scores of two or more 
sequences are equal, then a decision is made according to 
a more detailed comparison of the scores in SCORE (I). 

d. Subsequent Steps 

Since the most recent paths and scores are now 
in SSEQ(J) and SCORE(J), the roles of present and next are 
swapped, (I)«— (J), and step 3 is carried out again. Thus, 
after each decoding bit decision, the (I) and (J) are 
swapped, and step 3 is repeated time and again. Figure 15 
shows the decoding steps for four more shifts of the decoder 
input. The tables, after the eighteenth shift, are shown 


En figure 16. 
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It is noticed, at this time, that the first 8 
(DCL + K-1) bits are zero. This is explained by the pre- 
sense of 6(DCL) zeroes in the initial SSEQ(I) and the 
2(K-1) zeroes in the encoder when the first nonzero code 
bits were generated. 

Another note should be made of the ease with which a 
soft decision demodulation scheme could be incorporated 
into the scoring process. This dimension would produce a 
more accurate representation of the likelihood of a sequence 
and possibly eliminate the need for arbitrary decisions, 
which may also be time-consuming. ⁄ 

At a first glance this procedure may seem strange or 
awkward, but' there can be no doubt that the Viterbi algo- 
rithm is easily computer (hardware or software) implement- 
able. The recursive nature of the primary decoding step 


is the prime factor in controlling the size of that imple- 


mentation. 
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III. COMPUTER PROGRAM 


The term simulation in the title of this thesis should 
not be construed to mean the encoder and decoder operations 
(Programs) are simulated. By no means is this the case, 
for the rate 1/n convolutional encoder and the Viterbi 
decoder are software implementations that could be used 
in an actual system. Of course, some form of synchronization 
is needed for practical operation, but in this program that 
is assumed to have been accomplished by elements (hardware 
or software) preceding the decoder in the receiver. 

A. CHANNEL NOISE 

Actually, the simulation occurs when a channel is des- 
cribed in the program by a noise generating section. This 
section includes a quantization segment, random number seg- 
ment, noise generation segment, and perturbation (summing) 
segment. The following paragraphs are discussions of these 
four divisions of the channel noise program. 

l. Quantization 

This part of the program uses the encoded sequence 
as an input, then passes this sequence on to the noise sec- 
tions after Q zeroes or ones have been substituted for 
each 0 or 1. When Gel Les a program input this segment 
is bypassed for apparent reasons. 

2. Random Numbers 

The random numbers generated in this program were 


obtained by using the Lehmer congruential method. 
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Kr = aX, + b (mod Tg): 


where a = 257, b = 1, and To = p The term Xo, Starting 
number, is varied to provide representative sequences of the 
distribution shown in figure 18. These numbers are used to 
determine the time between set bits in noise sequence. The 
period of the random number sequence is 2.0. therefore, 
a large sampling of the sequence approaching 270 vould 
result in a binomial distribution. Hovever, the test runs 
used to determine results in this thesis employed message 
inputs of 10,000 bits in length. Thus, the distribution 
(dots) shovn in figure 17 is a better representation of the 
time betveen set bits in the noise sequence. 
3. Noise Generation 

After a random sequence has been generated in the 
previous segment, 2., a noise sequence is formed using 
the numbers. Multiplication of the number of set bits in 
16 bits (T), the length of one random number, with powers 
Of 2(2 ”, 2 1 20 əl, determine hov many zeroes vill occur 
before the next set bit in the noise sequence. For example, 
if the random number being considered by the program is 
0110011100100100, then its weight is 79 or 111,. If the 
noise parameters (core locations 10230 and 10236) specify 


multiply by PE then the number 7, is shifted one place 


8 
right and becomes er Likewise, if multiply by 2 is 
entered, then 79 becomes 165. These new numbers are then 


used to determine the space between set bits in the noise 


sequence. Each 16-bit number from the congruential random 
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number generator is used to generate a space. The number 
of 16-bit numbers used in one run is determined by the num- 
ber of code bits which are present when the noise is added 
to ES code sequence. 

The variation in the set bit density for the noise 
sequence is analogous to varying the density of channel 
errors. This fact is used in the evaluation of the codes 
in section IV. Figures 18 through 23 show the probability 
of k errors in 16 channel bits (Q = 1,2,3) versus the 
number of errors, k, in 16 channel bits. The plots are 
for averages (AT) of about 10 and 2 errors per 16 channel 
bits. Also shown is the corresponding plot of the Poisson 
distribution as a continuous curve. Any discrepancy between 
the theoretical and measured (dots) is attributed to the 16 
bit length constriction primarily. 

4. Perturbation 

The purpose of this last segment of the noise 
simulation is to simply exclusive OR (additive noise) the 
noise sequence to the quantized encoded sequence. This 
operation forms a perturbed sequence, which is the received 
sequence for the decoder. If this was to be changed to 
‘multiplication or some other operation another noise form 
could be simulated for the coding scheme. 

All of the details, as to programming specifics, are 


discussed in Appendix B, along with a program listing. 
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Figure 21. Simulated Noise (Q = 2, AT = 2.0) vs 
Poisson Distribution (AT = 2.0) 
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B. PROGRAM FLOW 


The four major sections of the program are the encoder, 


decoder, noise generator, and the message generator 


(ASCII 7-bit code). An appendix is devoted to a detailed 


listing and brief discussion of each of the above sections. 


The parameters that are needed to initiate the simula- 


tion are: 


ie 


N; the inverse rate (1/N +) (memory location 10200 


in Appendix F). 


K; the encoder constraint length (memory location 


10206 in Appendix FO. 

Q; the quantization levels (memory location 10214 in 
Appendix F). 

DCL; the decoder constraint length (memory location 
10222 in Appendix F). 

Noise parameters; operation (-/x) and operand 
(1,2,3,...) (memory locations 10236 and 10230 
respectively in Appendix F). 

Generator sequences; representations entered so 
that first cell in encoder is bit 0 (memory 
locations 13760-13776, where G, (13760, G,(13762,..., 


Gg (13776) ). 


Another programming concept, position independent code 


(PIC), 


is used to allow the user to move the entire program 


or any part (with few changes) to another part of core. 


This approach introduced the need for a program segment that 


initializes all addresses and counters (constants) used in 
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the program. The listing and a brief description of this 
segment are included in Appendix D. 

The entire program, in block form, is depicted in 
figure 24. If the reader is interested in a detailed 
program description he is referred to Appendices B-E. 
Appendix A contains a map of that part of core used for the 


implementation of the position independent code. 
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(10000-10170) 


Message Source 
ASC II / Bit Code 


Address and Counter 


Initialization 


(10710-11176) 


Rate 1/n Convolutional 
Encoder 





(11200-11706, 





Noise Simulator 


(11710-13176) 


Viterbi Decoder 


(13700-13400) |, 


Message Sink 


ASC II 7 Bit Code 





Note: The numbers in parentheses are 


actual addresses occupied by 
the instructions corresponding 


to a particular flow chart 
block. 


Figure 24. Block Diagram of Program Flow. 
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IV. RESULTS AND CONCLUSIONS 


The term best code is used in many papers, but for the 
most part the use of the term best code is not defined in | 
sufficient detail. This is one descrepancy that will be 
avoided by this thesis. Included in this section is a 
summary of the computer results obtained from a best code 
(defined below) determination procedure and encoder/decoder 


parameter variations. 


A. BEST CODE DETERMINATION 

The determination of a best code for a specific channel 
a simulated channel in this paper, follows a procedure based 
on Shannon's fundamental theorem for a discrete channel with 
noise. This procedure is best described using Shannon"s 
representation of the attainable region in a graph (figure 
25) of H(x) (information rate) versus Hy, (x) (probability 
of decoded messager error). If H(x) < C (channel capacity), 
then Shannon shows that 575 can be made arbitrarily small 
with a proper encoding procedure. Vhen H(x) » C, then the 
excess information being pushed onto the channel can only 
increase the uncertainty S of the decoded message. 
The minimum value of m is very close to H-C in thas 
latter case. 

The computer program in this thesis is implemented with 

a supplemental segment which stepped through all well-defined 


convolutional code generators of a specific code rate and 
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constraint length. Well defined refers to the constraint 
length definition of the code, meaning that the coefficient 
of the zero power in the generator polynomial must be set for 
at least one generator sequence, and then another sequence 
must have the coefficient set for the K-1 power. An 


example of this is shown here: 


n=2,K=4; g (D) = 1 + p? 
9, (О) лр. 
n= 3, K= 3; g (D) =D+ Des p 
g, (D) 2575-5055: 
93 (О) 2) oe ey DER 


"The rate is defined by the number of nonzero generator 
sequences. 

The first two paragraphs describe a procedure and the 
selection of encoders to be processed by the program. The 
program is applied to find the generator sequences, for a 
given code rate/constraint length, which come the closest to 
the lower error Ey си boundary of figure 25. Code rates 
of 1/2 and 1/3 are tested and if the input message is 
assumed to have 1 bit of information in each bit entering 
the encoder, then these rates correspond to those plotted 
on the horizontal scale of figure 25. Tables II through IV 
present a summary of the test run results for a best code 
determination. 

Each test run was made for a 40,000 bit encoder input 
applied over various samples of the simulated noise and at 
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Measure of average error 





ATTAINABLE 
REGION 


Tr formation Input Rate 


C = Channe1 Capacity 
Hu ax — Me 1mıized HİÜX) 
ae Transmitter 


Figure 25. The measure of average error for a given 
iNTfOTmatılmı input rate to a channel 
(figure is similar to Shannon"s 
Representation). 
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four values of AT(1,2,4,10 errors per 16 bits). The com- 
puter output is a listing of the decoded message errors for 
each well defined encoder. This output would require many 
additional pages to list in this thesis. Therefore, the 
summary in Tables II through IV is presented instead. The 
number of decoded errors is then divided by 40,000 to give 
the measure of average error for each test run (AT). 

The various values of AT are obtained by the noise 
simulation and have the characteristic distributions of those 
curves presented in figures 18 and 19. 

A discussion of these results is presented in part C of 


this section. 


B. DECODER PARAMETER VARIATIONS 

In order to provide an insight as to the effects of 
quantization (Q) and decoder constraint length (DCL), tables 
V and VI are shown. The effect of increasing Q or DCL is to 
decrease the frequency of errors thus improving the error 
correcting scheme. The codes used in these tables are those 


chosen as the best codes from tables II through IV. 


C. DISCUSSION AND CONCLUSIONS 

The determination of a best coding scheme for a specific 
channel application can be made from tables II-IV. The 
generator sequences which have the minimum frequency of 
errors (in place of Н (х)) for the various values of AT 
are outlined boldly. Using this comparison technique plus 


a comparison of the general statistics of the encoder at 


7] 





other values of AT, the following generator sequences are 
chosen to be the encoders which meet the criteria set forth 
previously. 


N= 2, K = 3; g, = (110), g, = (001) 


N 2, Ke 4, g, = (1110), g, = (0101) 
N= 3, K= 3; g, = (111), g, = (001), 
g4 - (001) 


These encoders were also used for the results obtained by 
varying Q and DCL in tables V and VI. 

With actual concrete measurements of the effectiveness 
(ninimum error) of a code and a means (computer program) of 
applying these codes to actual channel recordings there is 
no reason to have to settle for an inadequate error correcting 
system. The results presented in this section indicate clear 
preferences in choosing certain codes to accomplish desired 
communications in a given noise environment. The user speci- 
fics the task and with the computer a code can be chosen. 

The true value of such a program should be evaluated when 
implemented with actual channel noise. No one has been able 
to find a relation or algorithm for encoding that would enable 
a communication system to reach the maximum average informa- 
tion rate (C) with an arbitrarily small frequency of error. 
However, computer analysis offers the channel user an oppor- 
tunity to find a coding scheme which meets the requirements 


for that channel. 
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APPENDIX B 


PROGRAM LISTING OF NOISE SIMULATION 


The following flow chart depicts the flow of the 
instructions in the following machine language listing 
of the noise simulation subprogram. The numbers (base 8) 
on the upper left of the blocks in the flow chart corres- 
pond to those instruction addresses of the subprogram 


mentioned in that block. 
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(11200-11324) 
This segment quantizes the encoded 
block, that is, encoded l's become 
Q l's and encoded O's become Q O's. 
The result is then placed in block 
specified at (13604). 


(11330-11374) 
This segment generates the random 
number sequence by the Lehmer Con- 
gruential Relation. These numbers 
are stored in a block beginning at 
the address specified at (13010). 


(11376-11546) 
This segment uses the noise para- 
meters placed in locations (10230) 
and (10236) to change the value 

of the number of set bits in a 16 
bit random number, so that the 
density of errors (AT) in the noise 
is varied. The noise sequence is 
stored in a block beginning at an 
address in (13606). | 


111550-11646) i 

This segment determines the noise | 

value total (13666) and the number 

of errors (13670) for a given 

message, 1.e., a value 2 4 is an 
error ror Qz3. 


(11650-11700) y 

This segment adds the noise block 
to the quantized block and stores 
the result in a perturbed block 
beginning at an address in (13610) 


| END OF SUBPROGRAM | 
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The following computer printout is the noise program 
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APPENDIX C 


Program Listing of Rate 1/n 
Convolutional Encoder 


The folloving flov chart depicts the flov of the 
rate 1/n convolutional encoder implementation. The cor- 
responding machine language program follovs the flov 
chart and the instruction addressed are indicated on top 
of the respective flov chart block to vhich they corres- 


pond in the subprogram. 


(10710-10726) 


This segment sets up the addresses 
of the message block (RO), the En- 


coded block (Rl), and the stack 
(R6) 


(10730-10772) 

This segment uses the extended 
arithmetic element (AC-MQ) to 
shift in a new message bit to be 
encoded with the previous entries 
within the encoder constant length. 


(10774-11062) 
This segment determines the bits in 
AC(177302) that are to be encoded 
and steps through the generator 
Sequence block computing the n out- 
put code bits. 
(11064-11110) 
This segment contains loop indica- | 
tors for the number of bits left in 
MQ and the message word count. If 
the latter is not exhausted, a loop 


is made to put a new l6-bit word in 
MQ(177304). 


| END OF SUBPROGRAM | 
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The following computer printout is the encoder implementa- 
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APPENDIX D 
Program Listing of Viterbi Decoder 


The following flow chart depicts the flow of the 
Viterbi decoding algorithm implementation. The correspond- 
ing machine language program follows the flow chart and the 
instruction addresses are indicated on top of the respec- 
tive flow chart block to which they correspond in the 


subprogram. 
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(11710-12042) 
COMPUTE TRELLIS 


AND STORE IN TRELLIS 
BLOCK (13612). 





(12044-12146) i 
QUANTIZE TRELLIS TABLE 


ENTRIES FOR Q > 1 





(12150-12206) 





INITIALIZE SSEQ(I) (13614) AND 
SCORE (I) (13616) 









(02210123695) 


DECODE OLDEST BIT IN THE 
SSEQ(I) SEQUENCE WITH 
THE MINIMUM SCORE 





(12370-12546) v 


SHIFT INTO MQ(177304) 
N BITS FROM THE PERTURBED 
SEQUENCE (13610) 





(12550=12772) y 


COMPUTE A BETWEEN THE N-BIT INPUT 
AND EACH N-BIT BRANCH CODE. ADD 
THE BRANCH A TO IT CORRESPONDING 


PRESENT SCORE. DETERMINE THE 
BRANCH WITH THE MINIMUM SCORE 
ENTERING EACH NEXT STATE. 





(12774-13072) 


CHANGE EACH SSEQ 
CORRESPONDING TO THE BIT 
LOST FROM THE PRESENT STATE | 


(13074-13166) | 


-—— J 


(*) 
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0117-4 
OTT eb 
011748 
811742 
411744 
411746 
811750 
811752 
811754 
6811756 
81176 
011762 
011764 
811766. 
1172706 
011772 
011774 
1177€ 
12888 
a12a8082 
az128084 
GREG? 
a128168 
12612 
6412414 
ulu 
0120280 
И12022: 
41224 
412026 
012028 
И12022 
8128z4 
"12026 
B12044 
4012042 
иН12044 
612046 


tion just discussed. 


(MIESOO 


"082814 
"801672801 
e DDR 
POHLEP° 82 
¿84666 
2016204 
О 01208 
Н 50Н02 
cansada 
¿A16r14 
àÁü0124 
"885185 
pie Tu. 
г040205 
28016244 
cablirsoa 
гор 
2006285 
:005506 
885314 
NEIE 
d HD ZIL a 
PEHECHE 
"30598512 
#BB572 d 
"885314 
raaireil 
"1671 
"0801622 
f "85722 
"8821685 

0262802 

“11712 
/UBLA,7? 
Ban. na 
ce — НЛ2Н22 
cacha: 
A 
25: 
d DDDd Ze 
2::67 
“ütü 
"81642 
/¿BRRBR144C2 
POLE, ae 
"Ba1542 


86 


041245040 
012052 
412454 
012056 
6120868 
a41280682 
012464 
94012066 
012020 
Н012022 
и012074 
и012026 
a121806 
812102 

üflziüd 

812106 
6121108 
81211 
012114 
812116 
121268 
81212 2 
B12124 
812126 
4012150 
412132 
üizisd4 
üizis6 
ülzidü 
012142 
a12144 
12146 
12158 
и12152 
412154 
12156 
012160 
Sie ioe 
12164 
Wie Lee 
aie i 7 
и12172 
012174 
и12176 
812200 
И12202 
и0112204 
012206 


The following computer printout is the decoder implementa- 


"016781 
^"^0881e€^58 
/016702 
”0015584 
016282 
”081614 
20162804 
"8881614 
ЫН5Н12 
A/aage218 
¿103443 

21312 
ae dn 
casadas 
esüüğiiz 
2885212 
-. 


00572 eo 
7816703 
e881554 
"8895818 
"16784 
cM oed 
"ameziü 
POPP? dh!) 
e8e2218 
POF PSHE 
eaas?zh 
”885242 
eily (156 
2012200 
(8142 
ares ol 
fAALd ed 
/016282 
“481422 
culeros 
2 0814280 
2017321 
"a85ü22z 
“E15784 
mu = 2 
"8885821 
caer Ze 
01514 
arrg 





612210 
12212 
812214 
812216 
612226 
И12222 
812224 
412226 
a1zzs8 
8122522 
a122z4 
ai1z2ze 
412240 
a12242 
812244 
412246 
64122586 
a122^52 
812254 
12256 
012268 
a12262 
6412264 
Ии12266 
612270 
212 ә 
a12274 
И12226 
812568 
Н12202 
412304 
412306 
012210 
0123512 
G12 314 
812316 
a12z28 
ula 
127324 
812526 
81220 
0122 < 
812524 
Н125356 
4123486 
12342 
012244 


4125346 


2016?701 
/881468 
7160712 
/8891276 
/0167204 
d DPI Zä 
/016?6? 
/8014680 
/881418 
551276” 
d DDRD AP 
”“üldlatd 
/B12 zx 7 
/1272804 
Дие 
088 28 
/8081224 
/ 8526. 
/BB1462 
vübcüsi 
“712660 
/881444 
88013862 
cate2z6r 
fHALSE4 
/812262 
28801256 
¿BBSSES 
/B81346 
5 
¿acia 

isis 
/882766 
/ 885722 
“HAE? PE 
7166767 
“8142: 
“Hibisco 
“Te në 
PUL See 
euğtzi4 
7HSEP16 
/8081404 
/18688684 
PHOACEHE 
/BBB2408 
e HDD dén 
/8085214 


87 


c 
CA 


e 
812 
г 


non cn cn on 
PA AA Pa 


812266 
612571 
bizi? £ 
812 <: 4 
Ufo WE 
012466 
612462 
012484 
012406 
1241H 
a12412 
124114 
812416 
a12428 
812422 
412424 
612426 
a1z24z8 
И12422 
812454 
012426 
012448 
812442 
012444 
012446 
12458 
natos 
12454 
812456 
012460 
812462 
12464 
812466 
a1247? hs 
0124722 


/0855267 
“HALE? 
ceüüibiz 
/885724 
/0855262? 
/881256 
О 8188= 
/B88gaaag 
/080240 
¿BR 82408 
Peter oe 

/088820 
d DH ZA 
^ü0181682 
Z66702 
ehüisz4 
esse” 
fice see 
fHLEPHE 
/881276 
"ANC 
ISS 
fierst4 
CHOSEN 
/881212 
260614822 
far FGA? 

2 088414 
Zara er 
ZBRBB2 B 
“üüiire 
pulse. 
sir?:üz 
”881156 
/әӦӦ225? 
firrsid 
“8162 
ZEE 
г 77: üz 

cantares 
Ле тб 
¿aid 
2801166 
Cordan” 
/001258 
/н01816? 
2/ 881246 
/ 2026? 





6125168 
812512 
012514 
61216 
412524 
4l252z 
a12524 
B12526 
a125:8 
012552 
8125254 
HLESZE 
81254 
"ə 
01254 
Gig se 
412550 
E eM 
8129554 
412556 
012568 
PP SEE 
12564 


612602 
и12604 
и12606 
1261 
Н12612 
012614 
012616 
и12620 


и12624 
и12656 
1264 
иа12642 
12644 
B12646 


"081186 
/LB88818 
¿HERPES? 
PHHL226 
etüizsü 
PHEEZE? 
“üğizzü8 
à"ügi224 
üBgagaag? 
гв00248 
/166?6? 
с001206 
8121 
leE Gr 
"88984126846 
¿unicas 
A 
d DDpDH/ 
f/HHL268 
"811267 
"ügi18?6 
PALE? AE 
27:.1: 
PRLS or 

250202 
uü1044 
"812885 
fHFP 456" 
CORLEZE 
fOASHE? 
#BB1B26 
"8895167 
a"8801826 
28016?0Н5 
/201866 
Lé: KB ge 
fOHL114 
f7HOLaLS 
caderas 
dec ltaol18 
77 

ATTE 
/BG0774 
fHOHECE? 
canarrá 


E 


POPP E1e 


88 


812650 
[A 
412654 
412656 
иН12660 
Н12662 
812664 
812666 
012670 
Н12672 
tizt”4 
4Н012676 
Hier Ge 
B12 7 2 
4012204 
alero 
a1c?18 
12712 
812714 
812716 
012720 
De m 
a12024 
н12726 


012 2246 
и122080 
0812252 
012254 
5 
4121 


с о 
pa pa 


To ə nə 


01 
a1 


ми 


» n» n 
-J -I A ^J 74d 74 E 


=) =.) =f oti oa © 


D 
— 
də ə 


127 
812776 
412060 
013002 
01300d 
Hi1zBHc&à 


e" 886248 
¿HEEE? 
"888756 
fHOLE66 
"818267? 
"88874€ 
CHEEPE? 
ZBRBR184RF 
dDDD AP 
2::776£€ 
"Güg?z4 
Н Н 262 
"816786 
8887274 
fHLEPE? 
г "Ae í 202 
гво20Н226 
г012005 
fA? 456? 
"aga? 
Cu. 7 
"888718 
fHASLE? 
^Gg87168 
CIT OT 
agaro 
PALE EP 
BBB OE 
гои0624 
Dé 
d DDDE ZA 
1 'Radese 
"2:132 
"aSg6ez 


— 
Y C T 
T 
^ 
ín 
с m (H rr, £n 
"I SJ A m S N 


Pë, 
Ei 
=J 
J ün 
= 


e DDHGA/ 
“PHOEBE 
2626? 
“О Н8668 
“Hüsü 
d HD Zd at 
d 2162 aS 
¿BnR 742 








613614 
a1x881z 
13814 
136016 
613620 
6413422 
015024 
Н015026 
AL ZAG 
613432 
13834 
1336 
6130444 
4134042 
0115044 
0152046 
a1za858 
105€ 
64136434 
j ss 
1x86 
6130462 
(1435164 
13068 
61:<07 4 
0120722 
40150?4 
илә иә 


килсә 
Н8646 
CHLELE? 
”8080644 
PAGE? E? 
/akardd 
0888636 
2016785 
С008622 
2816786 
RAPLE 
2511016 
"0886216 
"885216 
2880418 
PHLE FAS 
POOP aa 
fPHLEP aE 
888676 
16715 
PHROSPE 
“11116 
PHHESLE 
"8167285 
¿ute él 


2 29 ZE 


— 


89 


Et EUH 
813182 
81x14 
ES 106 
081311 
15112 
13114 
asi UG 
13120 
ETSI 
13124 
1155166 
413130 
ERT, 
613134 
8151856 
6413140 
ülsidz 
13144 
0114€ 
413136 
Bos 
1:124 
8132156 
1:16f8 
75166 
012124 
E1156 


"AH TEE ELE 
О Н8658 
Seefe 
àÁU088646€ 
esügizfZe 
2655218“ 
eutazaü 
AEE a € 
PHAHAS4 
881212 
“ula” 
ed HDD a 
2226 
“üntded 
“1800406 
PHLEP AL 
PHRASES 
Elers 
^U3G04258 
ies)? 
“HHO? GG 
201867841 
28808442 
AORE TH z 
¿atado 
Albert? 
)Á860714 
JU ggg 





APPENDIX E 


Program Listing of Supplementary Subprograms 


1. Initialization subprogram (10174-10706) 
This subprogram enters the necessary addresses and 
constants into those locations noted in Appendix A by *. 


The computer listing below is the initialization subprogram. 


aigi?4 / 8167868 H10278 /612221 
a18081?6 £Süizr”e H10272 /B€2788 
HiGcaaó /GOBdg2 ataz?4 /88086868 
81282 ^/8t2z4?4 B10276 /H18821 
6810264 /Н12?67? 0а10300 /6662600 
010206 280888 010302 20053404 
ataz1a rabzd4ra B1G83Fn4 /H18021 
a10212 /aG12?67 ai18x6 "816727? 
6414214 “8388081 0102 16 ¿UBIICGÉE 
Н10216 /485464 Н10212 2122204 
иЫ10228 1276? 010214 Х816Р2? 
0810222 /88802Н0 0105 46 ZUERST 
16224 728024868 в18220 /17?z86 
theese ¿ULeror 7 2812267 
Н102520 ¿UBA 016224 /122204 
aida2sz "802456 818326 z/a8x4z8 
Bites + (PLETE? и10330 /063200 
010256 ИНННВ Н 618732 /177784 
010240 /005452 A1AZZA 28810821 
010242 78810781 ataxsse /881673 

818244 HELTAL 4163 7 2 02342 
810246 /605514 LATAS SUPP TAS 
aij c2z58 /84180167? 818244 PALEPE? 
010252 /88525252 0106346 4177284 
uludas “Yee? hd G46 356 703346 
Glaze /"BR8,8048£ 818152 /EÉEYT7BB 
aigzeaü ^/8410167? 4104354 1773804 
01062 /,002244 156 718821 
a10264 "625868 atga8seg "R63780 
810266 /082482 010262 /177204 


90 





416364 
a16036€ 
8180x798 
n18z72 
8180-74 
G610 376 
410406 
41402 
18484 
a18046806 
616410 
1841 
8186414 
8141€ 
a18428 
414422 
410424 
610426 
1040 
410432 
4108434 
414476 
18448 
a1688442 
818444 
18446 
168456 
a18045z 
6104454 
6160456 
410460 
818462 
aladed 
1804€ € 
818470 
4104472 
14724 
016047 € 
6105406 
10852 
414554 
(145485 
aigsi 
GTI Tz 
8160514 
1851€ 
816852ü 
116504522 
8168524 
о1052 

0104534 
410532 
6414534 


2010021 
290652080 
ld CHA 
/810021 
га 25002 
2 805202 
2016?05 
200222?6 
"a8gesüz 
“HPP she 
2010262 
f/HUSS16 
200682 
6200 
2/010021 
”006202 
cunas 
fO1b26F 
fHASe?4 
2005202 
2012205 
го6 8002 
casa Ab 
fHLb6E4 
(HPC EWE 
d HDD 
eüiüzer” 
400:2254 
causas 
a"at826? 
Á"amzes 
2 65202 
гин6202 
“810267 
2 080:2242 
Ди др 
CHREAG 
EE echte 
zül sy 
Z BB < 1 7 6. 
I C6 
fide E. 
eirrxüd 
2002206 
"BEP s f 
/42053164 
/107204 
7: 
‘Maxis 
/1722 e 
U13767 
Fa cs itd 
Ы Н8216Н 


91 


0180536 
610546 
8180542 
8180544 
818546 
616556 


016552 


610554 
818556 
a1856€8 
HLO562 
810564 
1856€ 
8185?8 
ИЛ0572 
61057 

A105 °6 
@16664 
010602 
ülücüd 
8168686 
6106164 
4108612 
18614 
1861€ 
818628 
И10622 
и10624 
ülüt 26 
aigez 
018652 
41084 
0180656 
410644 
414642 
610644 
414846 
n168656 
a18652 
1865354 
1656 
010666 
61066 
1686866864 
818666. 
a180678 
ALGER SY 
180674 
Ga10676 
ALA’? aA 
ALP be 

10284 
ninas 


/8016202 
2003142 
CHE 
005202 
/406502 
vü. 
“Bas as 
ZB1BZ67 
2002156 
‘6452023 
20062052 
/818262 
cdH CID 
eülit”82 
“ützi8t 
7: 
"Bee? 
2002 164 
01026т7 
MUST 
dDIGdDa 
20802074 
A 
aÁgasgas 
PHBRSAS 
fHHEZES 
uru 
2810226? 
ru 
Hiera 
‚suzred 
2166202 
202256 
“ulüz6” 
“0210 
*016202 
2862752 
ich 
eüadez?4c 
culazo? 
eüğsüet 
PALHP aa 
6220 
IAS 
“010062 
20 825804 
Z BB 8 2 6 ? 
¿B 20842 
fAHE 26? 
2002016 
d DG 
)C80Bgsad4z 
f/HO00068 





2. Message Input Subprogram (10000-10144) 

This subprogram was used to indicate that an alphanu- 
meric symbol typed at the Keyboard was entered into core 
(echo). The symbols are stored in core in an ASC II 
(7-bit) code representation. The Keyboard symbol, @ , 
is used to terminate message entry and the number of 
8-bit computer bytes used is stored in location (13710 
The following computer printout is a listing 


Appendix A). 


of the ASC II message input subprogram. 


818888 ^/818768068 ~ 8618864 /680212 
81G0882 AE И10066 /105727 
01860804 /88422 ülüt ”ü 7177564 
8168686 EE 160? 72180875 
8108810 /1B57 37 и010024 711412737 
018812 /172568 ül "/8880288 
18014 ^18z?25 01641080 Z477?566 
GUN “is LO 8418041602 /Zü??407 
GE /1:7 9686 atatad /425222 
Bibaz2 ¿122r10 0101026 27177564 
16024 ^/808z88 8181106 /T88375 
DEE E uf CIELO 
ültüsü /185722 01601414 Ber 12 
18802 /17 7564 aliie ¿Le 366 
818804 /16z?5 181268 EE 3 
MHİHMHZ6 7112037 81122 2812782 
HiDRAD /1?^?566 8160424 "868404 
a1884z2 /8g5282 416126 “1050206 
ülütdd 712372? иЛ0126 /065202 
atado 177562 Mais 7G’ rss 
илин Z/mün215 4106134 /806282 
04880592 /081:9256& 1816 "8868282 
416034 /112240 616140 "8168267 
a1ag5e /Za808248 18142 /£Büzü44 
818068 /7/1805?28 18444 OOOO 
TEE: (Ate Pot x: 





3. Message Output Subprogram (13200-13364) 

This subprogram is basically the same as that in 
section 2 757- appendix. Hovever, besides typing vhat 
is in core starting at a location specified at (13602), 
the number of errors (bit differences) are determined 
betveen the input message and the decoded out message. 


This message error value is stored at (13672) (Appendix A). 


The folloving computer printout is a listing of the ASC II 


(7-bit) code output subprogram. 


И12208 208167086 Gisger4 ¿106075 
8122082 “HEE Gee te dese 
6132604 ^;8012?01 ALS 3H 2868086812 
613206 / 0680048 013382 /17r56€ 
55-10 iiiü? G12<84 /105222 
Gisse Clg wee 13x86 CLI red 
012214 /122?28 бтн 2 M 
613216 /^/agB8086 812212 /28067?z4 
015228 /201455 012214 иб 0246 
1222 В PIHOSSS И15516 /000248 
5 firrsed hi3328 z88167?68t6 
И112226 7183785 8123222 86254 
“5 ea8a385z801 atssed 78167681 
015222 1808366 012226 0: 
01:5254 /1222180 20112220 805022 

613236 “OHHH 4115522 /811863 
012240 /801262 a1sss4 /874211 
4015242 filers? Н12226 /Н127804 
015244 78886815 613348 “0608086268 
Н12246 /127566 613342 /0686211 
013250 Cote rhe 413344 /8055822 
4812252 080812 613346 z/ür7?74&8z 
Ки Pha? Se Ra1zz58 ,/üazz8z21 
И12256 /17?564 4н12252 782868687 
QS 00 RES d 012254 "868822 

izz CLASE 137356 / 881365 
HİS “177:4:36£E 413360 /010262 
ALS266 /822206 a1xxeé2 / 06266 
Пи. uH Litas, 613764 dDDpppp 
ИН122722 /122564 " 


ә 





4. Analysis Subprograms (13250-13544) 

From address (13250) to address (13376) a short sub- 
program is listed that is used to step through well defined 
encoders of rate 1/2, 1/3, and 1/4. In the last addresses 
(13400-13544) a program is listed which determines the 
number of error bits in 16-quantized bits and stores the 
results so that a distribution of errors may be plotted 


as in figures 17-23. (See next page.) 
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41<250 
a13282 
a1z254 
815256 
Н1226Н 
015262 
012264 
412266 
Н12270 
ИН122722 


4Н12504 
0812206 
1535108 
ia dee 
612314 
b13316 
012220 
4812 
41 
T 
013334 
pm 
112332 

5. 
b13346 
üissaz 
613344 
61:46 
613356 
11:56 
a1zz&hd4 


[du L.J 
L^ m. 


cee 
ed 
226 


0153780 
UH xr 
ülssr?d 


maaş 


0133776 


/88526? 
Z/088512 
⁄ esL 
fHHHSHE 
fHHOH14G 
HLHH 
7812 767 
¿BRB88084 
7808088474 
ов 5267 
PHOAdES 


/002462: 
eadgaa018 
808168017 
fHL1E° 6? 
ÁAagaaad 
2/008458 
PHASE? 
/000442 
ALE er? 
/000426 
2002 
àÁaUa01606806 
zü” 
Á)oO08E512 & 
2008424 
2 862767 
/uhnadz 
2 80080414 
ими? 
dDDD Za ai 
üÜg 80416 
fHeer es 
ВВВ 02 
/handid 
сист 
4НА2?4 
0621 
*801402: 
ec DDDI Ze 
e488088 
cDDHI Ze 
UL e 
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6134066 
0152402 
6134064 
813486 
613410 
и15412 
012414 
012416 
01542Н0 
012422 
a1x424 
B1Z426 
4Н15420 
a1*x432 
813434 
815426 
613440 
012442 
015444 
413446 
6413450 
6413452 
613454 
6413456 
613460 
Hisdec 
Hisded 
413466 
61:40 
015422 
b1<4rd 
4815426 
613506 
4815502 
413504 
641546 
41351 (6 
413312 
a1s514 
AE 
[413528 
013522 
81:924 
MİZ 
A 
иИ12522 
0125 
diz 
41 
alsa 
ET Gen 


Le PO 


” 
> . 


3 
2 
4 
4 
4 


t İZİ td 
Cn сл 
pa no mqi 


4 


a 


7@16766 
“üunüzüz 
28167841 
/0080226 
2/0806281 
2808248 
¿duleror 
)UBRamaa 
PHA 
”?eülzrüd 
202416060 
a 
e DDDD AP 
224 
2 8002480 
/080240 
/280248 
г2058086 
2/816282 
2/0002:4 
280858082 
2006202 
20 06518 
usaste 
OS 
¿avala 
AÉ00147 
(OPP? 247 
20280201 
/8002401 
/ЫН5206 
288053825 
7 881414 
205887660 
d 2 005720 
2 2620 
“ütü” 
PHOALAES 
üg 888586 
/008240 
/ 6824 
Ant. 
/ 0 0002 
CHAHINE 
/080256 
/0865056 
2060684 
085214 
/000240 
Н 22? 
¿cacao ac 





RFFENDIX F 


SAMFLE EUN 


IHE CØOMHFUTER LISTING RELON THIS PRRAGEGER 1S A 
DEMONSTRATION GF THE PROCEDURE REGUIREG TO USE THE 
ENTIRE FROGRAM FROFERL Y. THE ELUBEHU IS STORED CN 
A DISK UNDER THE NAME VSAC SAW. THE CONVOLUTIONAL 
COLE USEC IS A RATE dea. (k=). CODE, OTHER INPUTS 
ARE Hei, DCL= 346 (¿26 IN EASE €), AND THE GENERATOR 


SEQUENCES GF 3 AND 4 (ERSE €). 


SEF VSAC. SAY 


= inet 11702 
Ur vai-u1 
Ai 
*xigzBamBaasadgc 
HLA She PGs 
kiBZ1478088881i 
*1022228888286 
*102502008001 
*102268888281 


s 
+123346100064 
x: 


е12г622000085 
"P ww 
"E 

. STHET 188868 
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The paragraph below is the input message for the sample 


run. 


THIS 15 A SAMFLE FUN OF A THESIS COMPUTER 

PROGRAM USED TO EVALUATE THE PERFORMANCE CF FATE 

APN CONVOLUTIONAL CODES, OVER A SIMULATEC NOISY 
CHANNEL, WITH RH VITEREI DECGODEEF. İTE (uc EHh  ETEKS 
FÜR THIS FUN CCORE LOCATIONS 1623260 AND 10336) ARE 
LHENSEIMEEURHSCOMERESUSEL TA OETRİH THE DİSTETEULTUN 
IN FIGURE 2Z@ (F THE THESIS TEXT. IHIS MESE ELL 
NOW EE TERMINATED EY THE “AT” CHARACTER ON THE FANEL. 


When the message above is encoded (as given above) and the 


noise is added, the resulting decoded message is: 


TAI FE RUNAS E? THESIS 
COMPUTE PROGEAM USED Nü EVRLURTE 
THESEERGNEURNCESUF RATE 1-H CUHVULUTTUNRLİBLUDES, 
BITE RAF TEET EE NOISY Фон ә 
El EEC Les A A ELE 
A CREER IEEE n 
Hite Soi Sie RUHUN PL E LS EHT {I 
CETRIN THE @ISTRIEUTION INEFIURE 
cL IHENLHECIPEENESD: (MiS MESSAGE 
HELL ə TEPE TİHETEU EL TEHH R" 
CHERCTER ÜN THM GENEG. 


When no coding is applied to the input message and the 


same noise is added the received message is: 


BLREGHIHBEKGBEERJSGMESVUN"HDROCTICEUCRCOMNU 
OR EORMCURMBED СИТСИЈАЈЕ DCGRBBDPEFÜERFED"F 
SKUI"R!REHVFVEH"EGCOERZGGNGD! NNHSPU"TCKUOI 
ULIGESONKY 15: 
FIIS zUG'HCORESLEITIOGJL A48 :?"Е0"942)4Ҹ95 
KaeEERIOUZIE:GEMGECUDGTOGCHSVEIN I DUDETFHCHUI 
VOREBZESGTRUISOLESHNOG VE TÜR 
ЕЈЕЈКЕРАК" NOUMIANEL. 8 f 
* Improper carriage return received. 


Zb 
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